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Abstract 

In  Argentinean  Patagonia,  there  exists  a  growing  interest  in  understanding  local  problems 
associated  with  the  increasing  of  wind  energy  penetration  in  the  grid.  Recent  papers  show  that  some 
of  these  power  quality  problems  can  be  reduced  by  the  proper  control  of  the  wind  farms.  In  this  way, 
this  work  deals  with  the  impact  and  potential  contribution  of  wind  generation  on  the  damping  of  the 
electromechanical  oscillations  called  inter-  and  intra-area  oscillations.  For  gaining  qualitative 
insights  and  understandings  on  this  complex  subject,  a  test  system  that  allows  to  isolate  the 
oscillations  modes  is  considered.  The  analysis  shows  that  new  control  concepts  for  wind  farm  can 
efficiently  contribute  to  the  power  system  damping. 
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1.  Introduction 

From  early  times  many  countries  have  been  modifying  their  legislation  in  order  to 
promote  the  development  of  renewable  energy.  Then,  political  and  economical  incentives 
tend  to  encourage  active  power  generation  as  much  as  possible  from  wind  farms.  As 
examples,  green  certificates,  different  kinds  of  tax  exemptions  and  financial  helps  to  the 
wind  power  generation  are  applied  in  Chubut  State  (Argentine)  where  the  biggest 
argentinean  wind  farm  is  located. 

However,  as  the  number  and  size  of  wind  farms  increase,  the  interaction  between  wind 
farms  and  grids  becomes  more  important  [1-4].  Then,  associated  with  the  total  amount  of 
power  connected  to  a  distribution  network,  some  power  quality  problems  can  arise  due  to 
the  fluctuating  nature  of  the  wind,  e.g.  voltage  and  power  disturbances  which  deteriorates 
the  network  performance.  This  fact,  which  has  limited  the  wind  farm  installations,  is 
potentiated  in  weak  grids  in  the  southern  tip  of  South  America  in  Patagonia,  where  the 
electrical  network  is  not  connected  to  the  main  interconnected  argentinean  power  system 
and  where  the  medium  wind  velocities  can  be  over  12  m/s  [5]. 

These  characteristics  have  generated  in  Patagonia  the  necessity  of  studying  the  impact  of 
wind  generation  on  electrical  networks  and  researching  about  the  potential  contribution  of 
wind  farms  to  overcome  some  network  problems.  In  this  way,  this  work  resumes  studies 
which  show  how  the  wind  farm  control  could  contribute  in  one  of  the  aspects  of  power 
quality:  the  damping  of  electromechanical  oscillations. 

Oscillations  are  undesirable  because  they  limit  power  transfers  on  transmission  lines 
and,  in  some  cases,  induce  stress  in  the  mechanical  shaft  of  synchronous  generators  [6]. 
These  oscillations  are  of  two  types  [6,7]: 


•  Local  mode  or  intra-area  oscillations.  They  are  associated  with  rotor  angle  oscillations, 
of  a  single  generator  or  a  single  plant  against  the  rest  of  the  power  system.  They  usually 
have  frequencies  from  1  to  3  Hz. 
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•  Global  mode  or  inter-area  oscillations.  They  involve  oscillations  of  a  group  of 
generators  in  one  area  against  a  group  of  generators  in  another  area.  Inter-area 
oscillations  are  in  the  range  of  less  than  1  Hz. 

In  conventional  networks  (i.e.  without  renewable  energy  resources  contribution),  small 
signal  analysis  and  control  are  well-established  techniques  [6-8].  These  techniques  are 
employed  in  this  work  looking  for  obtaining  indexes  which  allow  to  evaluate  the  potential 
contribution  of  wind  farms  for  damping  electromechanical  oscillations.  The  analysis 
considers  different  percentages  of  wind  power  penetrations  and  different  kinds  of  basic 
control  strategies.  Simulations  are  presented  to  support  the  analytical  results. 

The  structure  of  this  paper  is  as  follows.  In  Section  2,  a  wind  farm  model  is  presented.  In 
Section  3,  the  wind  farm  model  is  included  in  the  reference  frame  of  a  power  system  model. 
Section  4  presents  the  wind  farm  control  which  considers  different  strategies,  particularly  it 
includes  simple  control  strategies  to  regulate  frequency  and  voltage.  In  Section  5,  the 
impact  of  wind  farms  in  a  power  system  is  evaluated  through  an  eigenvalue  study,  and 
simulation  results  are  carried  out.  Finally,  in  Section  6  conclusions  are  presented. 

2.  Model  of  a  wind  farm 


A  complete  model  of  a  wind  farm  with  a  high  number  of  windmills  may  lead  to  compute 
an  excessive  and  impractical  number  of  equations.  The  size  of  the  wind  farm  model  may  be 
reduced  by  aggregating  several  wind  turbines  with  similar  incoming  wind  into  a  bigger 
turbine  called  aggregated  turbine  [9].  The  mechanical  and  electric  parameters  per  unit  are 
preserved,  and  the  nominal  power  is  increased  up  to  the  sum  of  the  nominal  power  of  the 
whole  set  of  turbines  to  obtain  the  parameters  of  the  aggregated  turbine. 

Making  a  decision  about  how  many  turbines  should  be  aggregated  depends  on  each 
particular  case.  Analyzing  network  changes  far  from  the  wind  farm,  it  could  be  enough  to  study 
the  wind  farm  as  an  only  one  aggregated  turbine.  However,  when  a  disturbance  takes  place  in 
the  wind  farm  it  could  be  necessary  the  study  of  all  the  generators  in  a  separate  way  [10]. 

As  the  main  objective  of  this  paper  is  to  examine  the  effects  that  the  wind  farm  has  over 
the  network  oscillations,  a  single  lumped  mass  of  the  aggregated  wind  turbine  is  considered 
[11-13]: 


h 

N2 


NQt  = 


J  &  + 


which  is  referred  to  generator  shaft  by  the  gearbox  transmission  N ,  being  Jt  and  /g  the 
inertias  of  the  aggregated  wind  turbine  and  the  generator  plus  the  gearbox  one,  and  Qt  and 
Qy  generator  and  turbine  speeds,  respectively.  Tt  is  the  turbine  torque  (Appendix  A): 


npr2 


v3Cp().)/Qt, 


and  Tq  is  the  generator  torque  defined  in  Appendix  B.  Because  of  in  Patagonia  most  of  the 
wind  generators  are  powered  by  double  fed  induction  (wound  rotor)  generators  (DFIGs), 
in  this  work  these  electric  machines  are  the  only  ones  considered,  even  when  the 
conclusions  can  be  applied  to  other  variable  speed  generators.  Due  to  the  electric  dynamics 
is  fast  against  that  one  of  the  rest  of  the  system  [14,15],  the  dynamical  aggregated  model  of 
a  wind  farm  can  be  approximated  by  the  corresponding  mechanical  behavior. 
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3.  Power  system  model 

The  linearized  differential-algebraic  (DAE)  model  for  a  conventional  multimachine 
power  system  (without  wind  farms)  is  [6] 


Ax  =  A  i  Ax  T  B\  A/g  T  P2  A  Eg  T  E\  Aw, 

(3) 

0  —  CiAx  T  P>iA/g  T  P^A  Eg, 

(4) 

0  =  CSAx  +  PE  A/g  +  D4A  Eg  +  PE  A  Ei, 

(5) 

0  =  P>6AEg  +  P>7AEi, 

(6) 

where  the  expressions  (3)-(6)  represent  synchronous  generator  dynamics,  static  stator 
synchronous  generator  equations,  active  and  reactive  power  equations  at  generator  buses 

rr 

and  at  load  ones,  respectively;  APg  =  [Ign  Iqn\  are  generator  currents;  A  Eg  =  [0n  Vn\ 

'T 

are  phases  and  amplitudes  of  voltages  at  generator  buses;  A u  =  [cos  T Mi  Erefi]  are  the 

T' 

synchronous  frequency,  mechanical  torques  and  reference  voltages,  and  AEi  =  [0n  Vn] 
are  phases  and  amplitudes  of  voltages  at  nongenerator  buses. 

In  order  to  include  the  wind  farm  in  the  electrical  network  equations  (3)-(6),  the 
linearized  model  of  the  aggregated  DFIG  wind  farm  operated  at  constant  wind  speed  is 
calculated  from  expression  (1): 


N?  ^  ^ r 


Ar‘urb  -  ATe  =►  co, 


N 


CD 


rnom 


( Ktwb/N2)(2/P)Acor  AT e 

■Pbase/cu  rnom  (2/P)  Pbase/&>  rnom  (2/P) 


(7) 


Eq.  (7)  is  expressed  in  p.u.,  being  Ht+g  =  1/2  ((/g  +  (Jt/N2))(wmom(2/P))2)/PnhSE,  the 
electrical  generator  speed,  ojrn0m  the  rated  electrical  generator  speed  and  Prase  the  rated 
power  for  scaling  all  equations  in  per  unit.  Also,  the  gain  Kt urb  is  obtained  from  Appendix 
A,  taking  the  derivative  of  Eq.  (16): 

_(npr2/2)v3ACp/AQt  -  (npr2/2)v3Cp] 

^turb  ?  V*/ 

Qii 


where  the  subscript  “1”  indicates  the  linearization  point. 

In  order  to  model  the  complete  system  (with  wind  farms),  it  is  necessary  to  modify 
expressions  (5)  (or  (6)  if  the  wind  farm  is  connected  at  a  nongenerator  bus)  with  the 
generated  active  power  and  to  modify  expression  (3)  as  a  consequence  of  the  new  state,  the 
wind  farm  equivalent  speed  (<x>r).  Also,  the  generator  torque  takes  different  expressions 
depending  on  the  control  strategy  applied.  Then,  before  proceeding  with  the  model  (APe), 
it  is  important  to  analyze  the  possibilities  about  the  control  of  the  active  and  reactive 
powers  of  the  wind  farm. 
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4.  DFIG  wind  farm  control  strategies 

As  it  is  pointed  in  Appendix  B,  vector  control  can  drive  independently  active  and 
reactive  powers  of  DFIG.  Then,  this  work  considers  the  next  cases: 

•  CP:  Active  power  control  (tracking  or  regulation)  and  reactive  power  control  AQe  =  0. 

•  FC :  Active  power  control  APe(A/)  (in  order  to  contribute  to  frequency  regulation)  and 
reactive  power  control  AQq  =  0. 

•  FCVC :  Active  power  control  APe(A/)  (in  order  to  contribute  to  frequency  regulation) 
and  reactive  power  control  A Qe(V)  (voltage  control). 


4.1.  Conventional  wind  farm  strategies  ( CP) 


One  of  the  most  conventional  modes  of  operation  of  wind  farms  powered  with  DFIGs 
has  been  to  track  the  curve  Tt0  (expression  (20),  Appendix  A): 


npr2 

~Y~ 


C  Q1 


thus  getting  the  maximum  active  power  available  in  steady  state:  Pi0  =  Ftof2t. 

In  this  operation  mode,  wind  farms  disturb  the  network  as  the  wind  velocity  changes. 
Then,  given  the  growing  contribution  of  wind  generation  in  electrical  systems  the 
aforementioned  strategy  of  maximizing  the  power  conversion  can  be  unsuitable  for  some 
cases.  Hence,  other  strategies  are  considered  such  as  the  wind  farm  power  regulation  at  a 
power  reference  (Psc)  [16]. 

Obviously,  no  contribution  to  reject  frequency  disturbances  can  be  expected  from  any 
two  previous  wind  farm  control  strategies.  This  is  because  of  the  fast  power  (vector) 
control  which  keeps  the  wind  farm  power  at  its  reference  value. 


4.2.  DFIG  wind  farms  with  frequency  control  (FC) 

In  order  to  develop  a  corrective  action  to  frequency  changes,  the  next  simple 
modification  of  the  power  reference  in  the  DFIG  wind  farm  is  evaluated  [17]: 

TV  =  Fsc  +  APe(A/),  (10) 

APm)  =  KAf\  (11) 

being  Pw f  the  wind  farm  power  reference,  P$c  the  static  power  reference,  and  APe^  a 
variable  component  of  the  wind  power  reference  which  is  proportional  to  the  frequency 
variation  A/. 

In  order  to  assure  the  proposed  control  strategy,  it  is  necessary  to  consider  an  inertia 
reserve  of  the  wind  farm  from  operating  each  wind  turbine  with  a  power  coefficient  Cpi 
lesser  than  the  optimum  one  (Appendix  A).  This  fact  implies  that  a  static  power  reference 
Psc  <i\o  must  be  assigned  at  the  nominal  frequency  /G.  Consequently,  the  wind  generator 
has  100  *  (Cpo  —  Cpi)/Cpo%  of  additional  power  to  contribute  to  the  frequency 
regulation  in  steady  state.  However,  it  its  important  to  note  that  during  a  transient,  the 
wind  farm  can  contribute  with  the  excess  of  energy  available  in  the  rotor  inertia. 
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Fig.  1.  Power-speed-frequency  curves. 


Fig.  1  shows  power-turbine  speed-frequency  curves  with  the  wind  velocity  as  a 
parameter.  In  particular,  wind  turbine  characteristic  curves  (expression  (17)  in  Appendix  A) 
are  depicted  in  the  power-speed  plane,  and  frequency  characteristics  are  presented  in  the 
power-frequency  plane.  Note  that  power-frequency  characteristics  are  equivalent  to  the 
speed  droop  characteristics  in  synchronous  generators  [8,7],  however,  the  slope  of 
power-frequency  characteristics  of  the  FC  wind  farm  can  be  more  pronounced.  In  this 
figure  the  operation  point  A  corresponds  to  the  P$c  for  nominal  frequency  fQ  and 
Cpi  <Cpo.  Because  of  the  turbine  speed  corresponding  to  the  point  A'  is  larger  than  the 
turbine  speed  in  point  B\  there  exists  an  excess  of  energy  stored  into  the  turbine  inertia.  On 
the  other  hand,  on  each  power-frequency  characteristic,  exist  two  upper  limits,  the  static 
limit  B"  and  the  dynamic  limit  C".  The  dynamic  limit  C ",  in  correspondence  with  Pto,  is 
specified  for  preserving  the  wind  farm  operation  when  the  frequency  remains  below  its 
nominal  value.  Then,  it  is  possible  to  use  the  energy  in  between  two  curves,  the  first  curve  is 
the  turbine  speed  characteristic  Pt  and  the  second  one  is  due  to  Pto  which  correspond  to 
Eqs.  (17)  and  (19)  in  Appendix  A,  respectively. 

In  the  present  work,  the  lower  limit  considered  for  the  generated  power  is  zero.  Then,  if 
frequency  increases,  the  wind  farm  can  reduce  its  active  power  injection  into  the  network 
allowing  that  the  energy  can  be  stored  in  the  inertia,  then,  all  the  excess  of  energy  goes  to 
the  inertia  reserve.  Hence: 

(12) 

Contrary  to  the  conventional  synchronous  generation,  double  fed  induction  machines 
are  able  to  modify  very  quickly  the  generated  power  in  comparison  with  the  rest 
of  the  electrical  system  [15].  The  aforementioned  implies  that  the  participation  of 
the  wind  farm  could  be  effective,  even  when  the  wind  power  involved  would  be 
significantly  lesser  (for  reasons  of  power  quality  [10])  than  the  conventional  installed 
power. 
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When  a  DFIG  wind  farm  is  operated  with  FC,  expression  (7)  is  modified  with: 

2 Ht+g  _  (-^turb/ N2)(2/ P)A<j)r  +  Kt^ Acor  . 

t^rnom  -^BASe/^  rnom  (2/P)  Pbase/<^  rnom  (2/P)  ’ 

being  A7,„t  =  -Psc/w^0  and  Kt01t  the  FC  gain. 

Then,  for  power  analysis,  expression  (5)  (or  (6))  is  modified  adding  the  linearized  expres¬ 
sion  for  the  FC  wind  farm: 

APe(A /)  =  KP(0'Acoc,  (14) 

considering  KPa)e  =  Kt0)Qm. 

4.3.  DFIG  wind  farms  with  frequency  and  voltage  control  (FCVC) 

About  the  control  of  the  wind  farm  reactive  power  A  ge,  strategies  accepted  nowadays 
are: 

•  Qe  regulation,  i.e.  A Qe  =  0  (assumed  in  previous  cases); 

•  cos  cp  regulation,  i.e.  cos  cp  =  specified  [18]. 

However,  taking  into  consideration  the  increasing  penetration  of  wind  farm  generation 
and  looking  for  satisfying  the  following  objectives  [7,19]: 

o  keeping  voltages  at  the  terminals  of  the  equipment  within  acceptable  limits, 
o  enhancing  system  stability  to  maximize  utilization  of  the  transmission  system, 
o  minimizing  the  reactive  flow  so  as  to  reduce  losses  ensuring  the  transmission  system 
operates  mainly  for  active  power  transfer, 

•  a  DFIG  wind  farm  voltage  regulation  is  evaluated,  i.e.  AQq  =  AQe( Vy 

Considering  active  compensation,  devices  as  synchronous  compensators  and  SVCs 
automatically  adjust  the  reactive  power  absorbed  to  maintain  voltages  of  the  bus.  From 
this  viewpoint,  DFIGs  with  reactive  power  control  can  also  be  considered  as  active 
compensation  devices. 

Taking  into  account  that  the  wind  farm  reactive  power  can  be  controlled  independently 
of  the  active  one,  expression  (5)  (or  (6))  is  modified  to  take  into  account  the  voltage 
control.  Here,  the  most  simple  case  is  considered: 

AQe(v)  =  KQVAV,  (15) 

being  KqV  the  control  gain. 

Note  that,  as  with  synchronous  compensators,  the  reactive  power  injection  will  be 
limited  by  the  nominal  apparent  power  Sn  of  the  DFIG  [18]. 

5.  Impact  and  potential  contribution  of  a  wind  farm  in  a  power  system 

5.1.  Selected  test  system 

In  this  section,  it  is  evaluated,  how  basic  control  strategies  of  the  DFIG  can  affect  the 
dynamical  behavior  of  a  power  system.  In  this  way,  Fig.  2  depicts  the  selected  test  system. 
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Fig.  2.  System  under  test:  inter-  and  intra-area  power  flows.  Thick  line:  intra-area  oscillation.  Dashed  line:  inter¬ 
area  oscillation. 


It  is  formed  by  two  areas.  Area  1  is  an  infinite  bus  meanwhile  Area  2  is  formed  by  two 
synchronous  machines  in  parallel  with  a  wind  farm  feeding  a  load  which  have  been 
considered  of  constant  impedance  type.  Each  synchronous  machine,  where  it  is  considered 
a  constant  motor  torque,  is  associated  with  a  standard  IEEE  Type  I  voltage  regulator  [7]. 

The  symmetrical  network  topology  considered  in  Fig.  2  is  specially  chosen  for  analyzing  in 
a  separate  way  intra-  and  inter-area  oscillations  [20].  This  topology  permits  to  simplify  the 
analysis  by  separating  the  calculus  of  the  network  eigenvalues  for  each  oscillation  mode. 

In  Fig.  2,  it  is  considered  that  the  disturbances  appear  in  the  loads  at  buses  1  and/or  Y. 
It  is  well  known  that  small  signal  analysis  allows  splitting  any  disturbance  input  in 
different  components  and  then  to  calculate  the  system  response  applying  the  superposition 
principle  [21].  To  the  ends  of  the  present  study  it  is  interesting  to  split  the  disturbance  input  in 
such  a  way  that: 


•  one  of  this  new  components  (named  signal  or  disturbance  or  fault  of  common  mode 
[22])  only  excites  the  inter-area  oscillation  and 

•  the  other  component  (named  signal  or  disturbance  or  fault  of  differential  mode  [22]) 
only  excites  the  intra-area  oscillation. 


In  this  way,  and  taking  into  account  the  network  topology  of  Fig.  2,  any  fault  at  bus  1  and/ 
or  Y  (Eli  and  PL1/,  respectively)  will  be  divided  in  a  common  load  signal  [22]: 


Tlc 


Eli  +  Pyv 


•> 


2 
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which  is  applied  at  the  same  time  at  both  buses  1  and  F,  and  a  differential  load  signal:  [22] 


P  LD 


^L1  Py\f 
2 


? 


applied  at  the  same  time  in  both  buses  1  and  F  but  with  opposite  signs  in  order  to 
reconstruct  Py\  and  PLy: 


P  Ll  =  P  LC  +  P  LD? 

Plv  =  Plc  ~  PhD- 

Note  that,  in  Fig.  2,  inter-area  oscillations,  and  no  intra-area  ones,  are  obtained  from  Plc 
(i.e.  the  component  of  fault  which  is  equivalent  to  two  equal  signals  with  the  same  sign  at 
buses  1  and  F).  In  this  case  both  generating  groups  (synchronous  generators  and  wind 
farms)  behave  in  the  same  way  respect  to  the  infinite  bus,  i.e.  the  synchronous  generator 
shafts  experience  the  same  motions  indicating  that  all  the  Area  2  behaves  coherently 
against  the  Area  1 .  Then,  the  power  flow  in  between  the  areas  corresponds  to  the  indicated 
in  Fig.  2  with  dashed  lines.  As  a  consequence  of  the  network  symmetry  and  the 
characteristic  of  the  Plc  component  of  the  disturbance,  there  is  no  flow  in  between  the 
synchronous  generators  (buses  1  and  F),  i.e.  there  are  not  intra-area  oscillations. 

On  the  other  hand,  intra-area  oscillations  can  be  studied  from  the  differential 
component  of  the  disturbance  which  perturbs  the  network  at  buses  1  and  F,  the  same 
quantity  and  at  the  same  time  but  with  opposite  signs.  Due  to  the  characteristic  of  Pld  and 
the  network  topology,  the  excess  of  power  in  one  generator  group  in  Area  2  is  exactly 
absorbed  by  the  lack  of  power  on  the  other  generator  group  in  the  same  area.  As  a 
consequence  of  the  power  flow,  which  is  depicted  in  thick  line  in  Fig.  2,  there  is  no 
variation  of  the  power  flow  generated/absorbed  by  the  infinite  bus,  i.e.  there  are  not  inter¬ 
area  oscillations. 

Finally,  the  real  response  of  the  system  is  obtained,  by  applying  the  superposition 
principle,  from  the  addition  of  the  system  responses  to  both  components  of  the 
disturbance. 


5.2.  Cases  of  study.  Simulation  results 

In  order  to  evaluate  the  power  system  dynamic  behavior  when  wind  power  penetration 
is  increased  (the  same  quantity  that  synchronous  generation  is  decreased),  three  cases  are 
studied  by  means  of  an  eigenvalue  analysis  based  on  Fig.  2: 

1.  wind  farm  with  constant  active  and  reactive  power  regulation  (CP),  i.e.  A PQ  =  0  and 
Age  =  0, 

2.  wind  farm  with  FC,  i.e.  A PQ  =  APe(A/)  and  A Qe  =  0, 

3.  wind  farm  with  frequency  and  voltage  control  (FCVC),  i.e.  A PQ  =  APe(A/)  and 

A  Qq  =  A  Qe(y)- 


Fig.  3  presents  simulation  results  for  the  CP  case  taking  into  consideration  the  selected 
system  in  Fig.  2  with  a  150  MW  load  at  both  buses  (1  and  F)  and  wind  farm  and 
synchronous  generations  of  30  and  120  MW  at  each  bus,  respectively.  A  disturbance  of 
40  MW  at  bus  1  is  considered  at  t  =  2  s.  The  transient  response  of  frequency  at  bus  1  is 
presented  in  part  (a),  this  response  is  separated  in  intra-  and  inter-area  oscillations  in  parts 
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a 


Fig.  3.  Dynamic  response  of  the  system  considering  the  CP  case  and  a  wind  farm  penetration  of  20%:  (a) 
frequency  at  bus  1,  (b)  intra-area  oscillation  and  (c)  inter-area  oscillation. 


(b)  and  (c),  respectively.  The  intra-area  oscillation  corresponds  to  the  difference  between 
shaft  velocities  of  synchronous  generators  at  buses  1  and  Y  and  the  inter-area  response  is 
obtained  by  taking  away  the  inter-area  frequency  of  part  (b)  from  the  transient  response  of 
frequency  of  part  (a). 

In  next  subsections,  the  effects  of  FC  and  FC  and  voltage  control  (FCVC)  cases  on  both 
kinds  of  oscillations,  intra-  and  inter-area  ones  are  analyzed  and,  finally,  the  total  response 
of  the  electromechanical  oscillation,  voltage  and  control  effort  are  presented. 

5.3.  Intra-area  oscillations 

Based  on  the  considerations  commented  in  Section  5.1,  the  eigenvalues  corresponding  to 
the  intra-area  oscillations  are  calculated  and  presented  in  Table  1.  Parts  (a)-(c)  of  Table  1 
corresponds  to  CP,  FC  and  FCVC  cases,  respectively.  Note  that,  looking  for  emphasizing 
the  problems  under  study,  the  eigenvalues,  which  are  scaled  by  2n  giving  frequency  in  Hz, 
are  calculated  for  a  broad  span  of  wind  penetrations. 

From  analyzing  Table  1  it  is  observed  that: 

•  At  the  same  wind  power  penetrations,  in  parts  (a)-(c),  the  most  significative  differences 
correspond  to  the  eigenvalues  in  between  lines. 
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Table  1 

Intra-area  oscillation 


Peol  0  MW 

Peol  30  MW  Peol  60  MW 

Peol  90  MW 

Peol  120  MW 

Peol  150  MW 

(a)  Case  1:  DFIG  wind  farm 

with  constant  power  and  Qe 

,  =  0 

-157.62 

-157.62 

-157.62 

-157.62 

-157.62 

-157.63 

-4.31 

-4.29 

-4.27 

-4.25 

-4.25 

-4.25 

-2.95 

-2.95 

-2.94 

-2.94 

-2.95 

-2.95 

-0.42 

-0.40 

-0.39 

-0.39 

-0.39 

-0.39 

±1.15i 

±1.15i 

±1.15i 

±1.14i 

±1.14i 

±1.141 

-0.075 

-0.093 

-0.108 

-0.117 

-0.122 

-0.123 

±1.42i 

±  1 .421 

±1.42i 

±1.42i 

±1.43i 

±  1 .421 

0 

-0.18 

-0.18 

-0.18 

-0.18 

-0.18 

(b)  Case  2:  DFIG  wind  farm 

with  frequency  control  and 

Qe  =  0 

-157.62 

-157.62 

-157.62 

-157.62 

-157.62 

-157.63 

-4.31 

-4.34 

-4.51 

-4.97 

-5.96 

-6.01 

-2.95 

-2.95 

-2.95 

-2.94 

-2.95 

-2.95 

-0.42 

-0.41 

-0.44 

-0.46 

-0.43 

-0.42 

±1.15i 

±1.14i 

±1.09i 

±1.211 

±1.181 

±1.171 

-0.075 

-0.17 

-0.38 

-1.05 

-1.05 

±1.42i 

±  1 .421 

±1.38i 

±1.49i 

±1.50i 

0 

-0.18 

-0.18 

-0.18 

-0.18 

-0.18 

(c)  Case  3:  DFIG  wind  farm 

with  frequency  and  voltage 

control 

-157.62 

-157.62 

-157.62 

-157.62 

-157.62 

-157.63 

-4.31 

-4.36 

-4.55 

-5.02 

-6.01 

-6.02 

-2.95 

-3.12 

-3.28 

-3.35 

-3.38 

-3.38 

-0.42 

-0.40 

-0.43 

-0.45 

-0.39 

-0.39 

±1.15i 

±1.13i 

±1.06i 

±1.211 

±1.05i 

±1.041 

-0.075 

-0.17 

-0.39 

-0.70 

-1.05 

-1.06 

±1.42i 

±  1 .421 

±1.37i 

±1.03i 

±1.49i 

±1.501 

0 

-0.18 

-0.18 

-0.18 

-0.18 

-0.18 

Eigenvalue  analysis. 


•  These  eigenvalues  present,  at  each  wind  power  penetration,  a  higher  damping  factor  for 
the  FC  and  FCVC  cases  with  respect  to  the  CP  one.  This  implies  a  potential 
improvement  in  the  damping  of  the  intra-area  oscillation  for  the  FC  and  FCVC  cases. 

•  For  each  control  strategy,  the  damping  of  these  eigenvalues  increases  with  the  wind 
penetration.  This  implies  that  the  damping  of  the  intra-area  oscillation  increases  with  wind 
penetration.  Again,  this  increment  is  much  more  marked  in  the  FC  and  FCVC  cases. 

•  In  the  CP  case,  the  eigenvalues  displacement  with  the  increment  of  the  wind  energy  is 
not  significative.  This  fact  is  attributable  to  the  independence  of  the  CP  strategy  with 
respect  to  the  frequency  and  voltage  variations.  In  this  case,  the  eigenvalues  movement 
is  a  consequence  of  unloading  the  synchronous  generators  with  wind  penetration,  i.e. 
there  is  a  drop  in  the  load  which  “sees”  each  synchronous  generator  provided  that  more 
wind  energy  is  injected  into  the  network. 
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Figs.  4  and  5  show  the  movement,  for  different  wind  energy  penetrations,  of  the 
aforementioned  eigenvalues  for  the  CP  and  FC-FCVC  cases,  respectively.  Fig.  5  indicates 
the  strong  (positive)  influence  on  the  damping  of  the  intra-area  oscillation  with  the  FC  and 
FCVC  control  strategies. 

Fig.  6  presents,  for  the  three  strategies,  the  intra-area  oscillation  when  takes  place  a  load 
change  in  40  MW  at  t  =  2  s  at  bus  1,  being  the  wind  farm  power  generation  and  the 
synchronous  power  generation  of  30  MW  and  120  MW  at  each  bus,  respectively. 
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Fig.  4.  Intra-area  eigenvalues  for  the  CP  case  at  different  wind  farm  penetrations. 
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Fig.  5.  Intra-area  eigenvalues  for  FC  and  FCVC  cases  at  different  wind  farm  penetrations. 
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Fig.  6.  Intra-area  oscillation  and  voltage  at  bus  1  with  a  wind  farm  penetration  of  20%.  Wind  farm  with:  constant 
power  (CP),  frequency  control  (FC)  and  frequency  and  voltage  control  (FCVC)  strategies. 


Simulation  results  agree  with  the  expected  improvements  from  Table  1  and  show  the 
strong  impact  on  the  intra-area  oscillation  of  the  FC  and  FCVC  control  strategies  when 
they  are  compared  with  the  CP  one.  At  the  same  time  it  is  observed  that  the  transient 
behavior  of  the  intra-area  mode  is  almost  the  same  for  the  FC  and  FCVC  control  strategies 
but  the  last  one,  as  a  consequence  of  the  reactive  power  loop  (Age(F)),  produces  a  better 
voltage  profile. 

5.4.  Inter-area  oscillations 

Table  2  presents,  based  on  the  considerations  commented  in  Section  5.1,  the  eigenvalues 
corresponding  to  the  inter-area  oscillations  for  the  three  aforementioned  cases. 

Analyzing  Table  2  the  next  conclusions  (similar  to  the  Table  1)  can  be  obtained: 

•  Comparing  the  inter-area  eigenvalues  for  the  three  control  strategies  at  equal  wind 
power  penetrations,  the  main  differences  correspond  to  the  eigenvalues  in  between  lines. 

•  Also,  for  each  case,  the  damping  of  these  eigenvalues  increases  with  the  wind 
penetration  and  it  is  much  more  marked  in  the  FC  and  FCVC  cases.  This  is  shown  in 
Figs.  7  and  8  where  are  presented  the  evolutions  of  the  aforementioned  eigenvalues  for 
the  CP  and  FC-FCVC  cases,  respectively.  This  implies  that  the  damping  of  the  inter¬ 
area  oscillation  increases  with  wind  penetration. 

•  These  eigenvalues  present,  at  each  wind  power  penetration,  a  damping  factor  which  is  highly 
superior  for  the  FC  and  FCVC  cases  with  respect  to  the  CP  one.  This  implies  a  substantial 
improvement  in  the  damping  of  the  inter-area  oscillation  for  the  FC  and  FCVC  cases. 
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Table  2 

Inter-area  oscillation 


Peol  0  MW 

Peol  30  MW  Peol  60  MW 

Peol  90  MW 

Peol  120  MW 

Peol  150  MW 

(a)  Case  1:  DFIG  wind  farm 

with  constant  power  and  Qe  - 

=  0 

-157.62 

-157.62 

-157.62 

-157.62 

-157.62 

-157.63 

-3.76 

-3.71 

-3.68 

-3.66 

-3.64 

-3.64 

-2.70 

-2.70 

-2.70 

-2.71 

-2.72 

-2.74 

-0.49 

-0.48 

-0.48 

-0.48 

-0.48 

-0.49 

d=1.46i 

dzl.47i 

dzl.48i 

dzl.49i 

dzl.50i 

dzl.50i 

-0.039 

-0.046 

-0.051 

-0.056 

-0.059 

-0.061 

dzl.OOi 

dzl.OOi 

dzl.OOi 

±0.99i 

±0.99i 

±0.99i 

0 

-0.18 

-0.18 

-0.18 

-0.18 

-0.18 

(b)  Case  2:  DFIG  wind  farm 

with  frequency  control  and  Q 

e=0 

-157.62 

-157.62 

-157.63 

-157.63 

-157.63 

-157.63 

-3.76 

-3.75 

-3.90 

-2.73 

-2.38 

-2.36 

-2.70 

-2.70 

-2.71 

-1.60 

-2.42 

-2.41 

-0.49 

-0.46 

-0.46 

-0.50 

-0.51 

-0.51 

d=1.46i 

dzl.48i 

dzl.53i 

dzl.53i 

dz 1 .5 li 

dzl.50i 

-0.039 

-0.20 

-0.58 

-4.49 

-6.27 

-6.38 

dzl.OOi 

±0.98i 

±0.74i 

&-1.60 

&-0.23 

&-0.23 

0 

-0.18 

-0.18 

-0.18 

-0.18 

0.18 

(c)  Case  3:  DFIG  wind  farm 

with  frequency  and  voltage  control 

-157.62 

-157.61 

-157.61 

-157.60 

-157.60 

-157.60 

-3.76 

-3.77 

-3.94 

-3.34 

-3.38 

-3.38 

-2.70 

-3.02 

-3.26 

-1.70 

-2.53 

-2.52 

-0.49 

-0.46 

-0.45 

-0.47 

-0.44 

-0.44 

±1.46i 

dzl.46i 

dzl.47i 

dzl.38i 

dzl.27i 

dzl.26i 

-0.039 

-0.20 

-0.58 

-4.51 

-6.26 

-6.37 

dzl.OOi 

±0.98i 

dz0.74i 

&-0.47 

&-0.23 

&-0.23 

0 

-0.18 

-0.18 

-0.18 

-0.18 

-0.18 

Eigenvalue  analysis. 


Simulation  results,  presented  in  Fig.  9,  show  the  inter-area  oscillation  and  the  voltage  at 
bus  1.  A  strong  impact  on  the  inter-area  oscillation  of  the  FC  and  FCVC  control  strategies 
when  they  are  compared  with  the  CP  one  is  observed.  Note  that,  the  inter-area  oscillation 
presents  a  slightly  better  behavior  for  the  FC  case  when  it  is  compared  with  the  FCVC 
case.  However,  the  last  one  improves  the  voltage  profile. 

5.5.  Complete  transient  response  at  bus  1 

Fig.  10  presents,  for  the  three  cases  and  at  bus  1,  in  part  (a)  the  complete  transient 
response  of  the  frequency,  i.e.  the  shaft  oscillation  of  the  synchronous  generator  at  bus  1, 
in  part  (b)  the  voltage  and  in  part  (c)  the  apparent  power  of  the  wind  farm. 
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Fig.  7.  Inter-area  eigenvalues  for  the  CP  case  at  different  wind  farm  penetrations. 


Fig.  8.  Inter-area  eigenvalues  for  FC  and  FCVC  cases  at  different  wind  farm  penetrations. 


Parts  (a)  and  (b)  of  Fig.  10,  clearly  show  the  improvements  over  the  transient  response  of 
frequency  and  the  voltage  for  the  FC  and  FCVC  cases  when  they  are  compared  with  the 
CP  case.  The  reactive  power  loop  (A ge(F))  produces  in  part  (a)  of  Fig.  10  a  slight  decrease 
in  the  performance  of  the  transient  response  of  the  frequency  with  respect  to  the  FC  case. 
This  fact  is  attributable  to  the  lesser  variation  in  the  voltage  observed  in  part  (b).  Indeed,  as 
the  voltage  changes,  constant  impedance  loads  (considered  in  this  work)  tends  to  damp  the 
frequency  oscillation  by  modifying  the  active  power  absorbed. 

In  part  (c)  of  Fig.  10,  it  is  noted  that  when  the  disturbance  appears,  DFIG  wind  farms 
with  FC  and  FCVC  react  carrying  out  a  contribution  of  active,  and  active  and  reactive 
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Fig.  9.  Inter-area  oscillation  and  voltage  at  bus  1  with  a  wind  farm  penetration  of  20%.  Wind  farm  with:  constant 
power  (CP),  frequency  control  (FC)  and  frequency  and  voltage  control  (FCVC)  strategies. 


powers,  respectively.  The  reactive  power  contribution  is  corroborated  in  part  (b)  where  is 
observed  the  improvement  in  the  voltage  profile.  Note,  that  the  maximum  (transient) 
increment  of  the  apparent  power,  i.e.  the  control  effort  involved,  is  about  a  fourth  part  of 
the  steady  state  condition  with  less  than  a  half  second  of  duration.  However,  this  is  enough 
to  strongly  improve  the  transient  response  of  frequency  in  part  (a)  and  the  voltage  response 
of  part  (b). 

6.  Conclusions 

In  Argentinean  Patagonia,  there  exists  an  increasing  interest  in  understanding 
some  problems  associated  with  wind  power  generation.  This  fact  is  motivated  for 
particular  characteristics  of  the  region  as  weak  grids  and  high  wind  speeds.  In  this  way,  this 
work  presents  studies  which  show  how  control  strategies  of  the  wind  farms  could 
contribute  in  one  of  the  aspects  of  power  quality:  the  damping  of  electromechanical 
oscillations.  Also,  considering  the  apparent  power  available,  a  reactive  power  control 
strategy  which  enables  to  improve  the  voltage  profile  at  the  bus  where  the  wind  farm  is 
connected  was  evaluated. 

Three  cases  were  considered,  wind  farms  with  conventional  power  control  (constant 
power:  CP),  with  frequency  control  (FC)  and  with  frequency  control  and  reactive  power 
control  (FCVC). 

The  analysis  of  the  impact  of  wind  power  on  the  electromechanical  oscillations  has  been 
done  by  observing  the  location  of  eigenvalues  for  different  wind  power  penetrations  in  a 
test  system  that  allows  to  isolate  the  oscillations  modes.  It  was  showed  that  as  the  wind 
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Fig.  10.  System  response  at  bus  1  with  a  wind  farm  penetration  of  20%:  (a)  frequency,  (b)  wind  farm  voltage  and 
(c)  apparent  power  of  the  wind  farm. 


power  penetration  was  increased  wind  farms  with  FC  and  with  FCVC  control  strategies 
highly  increased  the  eigenvalues  damping  for  both  kinds  of  oscillations,  inter-  and  intra¬ 
area  ones.  Then,  these  strategies  contribute  to  the  damping  of  the  electromechanical 
oscillations.  This  was  corroborated  with  simulation  results.  Also,  the  reactive  power 
control  allowed  to  regulate  the  voltage  with  a  slight  decrease  in  the  performance  of  the 
transient  response.  Also,  simulation  results  showed  that,  with  a  reasonable  control  effort 
(incremental  wind  farm  apparent  power),  is  possible  to  highly  damp  electromechanical 
oscillations. 

Finally,  note  that,  due  to  the  basic  characteristic  of  the  control  strategies  considered, 
it  can  be  expected  that  investigations  on  other  network  configurations  will  led  to 
similar  results,  counteracting  adverse  effects  of  high  wind  energy  penetration  in  power 
systems. 
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Appendix  A.  Wind  turbine 


The  torque  and  the  mechanical  power  developed  by  a  wind  turbine  are  given  by  [23] 
Tt  =  -^-v3Cp(X)/Qu  (16) 

=nf*ft3Cp(A),  (17) 

where  p  is  the  air  density,  r  is  the  radius  of  the  turbine,  v  is  the  wind  speed,  Cp(/ 1)  is  the 
power  coefficient,  X  —  Q{r/v  is  the  tip  speed  ratio  and  is  the  turbine  speed. 

Because  of  Cp  is  maximum  at  the  optimum  tip  speed  ratio  X  —  X0,  to  extract  the 
maximum  available  power  at  any  wind  velocity,  the  control  system  should  adjust  the 
turbine  speed  to  force: 

Qt  =  Qt0=—.  (18) 

r 

At  this  rotational  speed  the  maximum  turbine  power  Pt0  and  the  corresponding  torque  T t0 
result: 


npr2  ( Qtr\ 3 

2  UcJ  P° 


npr2 


CpoQ 


2 
t  9 


being  Cpo  =  Cp(20)  the  maximum  power  coefficient. 


(19) 

(20) 


Appendix  B.  Generating  units  (DFIG) 

Variable  speed  wind  farms  powered  by  DFIG  (wound  rotor)  are  the  only  ones 
considered  in  this  work.  The  rotor  winding  is  fed  using  a  back-to-back  voltage  source 
converter  as  depicted  in  Fig.  1 1 . 

The  operation  of  an  induction  machine  can  be  analyzed  using  the  classical  theory  of 
rotating  fields  and  the  well-known  d-q  model  [24].  In  this  way,  the  general  reference  frame 


Fig.  11.  Generating  unit:  double  fed  induction  generator  (DFIG)  with  rotor  converter  for  variable  speed 
operation. 
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model  for  the  induction  machine  is  expressed  through  the  following  matrix  equation: 


Usx 

Rs  +  pLs 

cogLs 

PLm 

ojgLm 

4y 

Us y 

cogLs 

Rs  +  pLs 

(jjgLm 

PP  m 

hx 

uYX 

PL  m 

(cOg  (ov)Lm 

Rr  +  pLY 

— (cog  —  ojr)Lr 

hy 

UYy 

((Ug  (jJy)Lm 

PL  m 

{o)g  —  ojr)Lr 

Rr  +  pLv 

_  irx  _ 

(21) 

where  Ls  and  Lr  are  the  stator  and  rotor  inductances,  Lm  is  the  magnetizing  inductance,  Rr 
and  Rs  are  the  stator  and  rotor  side  resistances,  cor  is  the  electrical  rotor  speed,  cog  is  the 
reference  frame  speed,  usx  and  usy  are  stator  voltages  and  urx  and  ury  are  rotor  voltages, 
finally  isx,  isy,  irx  and  iry  are  the  stator  and  rotor  currents.  If  cog  is  equal  to  the  line 
frequency  cae,  the  synchronously  rotating  reference  frame  model  is  obtained. 

The  electromagnetic  torque  Te,  the  stator  active  (Pe)  and  reactive  ( Qe )  powers  of  a 
DFIG  are  given  by 


Te 

Pq 

Qe 


9  PPm(jsxhy  4 y 4x) 


sr 


2^Usxhx  T  Usy  isy)  i 

9  (Usy 4x  Usxisy), 


(22) 

(23) 

(24) 


with  P  the  number  of  poles.  Because  of  vector  control,  it  is  possible  to  control  active  and 
reactive  powers  independently  [24,25]. 


Appendix  C 

Wind  turbine:  P  =  700  kW  at  v  =  11.75  m/s;  r  =  22m;  p  =  1.224;  Cp(2)  =  1.206  x  10-527 
-2.71  x  10-426  +  2.0399  x  1(T325  -  6.519  x  1(T324  +  1.257  x  1(T223  -  1.54  x  1(T222+ 
2.85  x  10-22  —  4.05  x  1(T2. 

DFIG :  V  =  690  V;  Rs  =  0.0067  Q;  P  =  4;  Ls  =  0.0075  Hy;  Lm  =  0.0061  Hy;  Lr  =  0.0084 
Hy;  Ry  =  0.187Hy;  a  =  0.3806  turns  ratio;  /t+m/  =  536096.59  Kg  m2;  N  =  51. 
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